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SECTION  I 


INTRODUCTION 

This  paper  presents  a  derivation  of  the  equations  of  motion 
of  a  gas  whose  molecules  may  be  modelled  as  rigid  rotors.  A 
rigid  rotor  has  three  translational  degrees  of  freedom  and  two 
rotational  degrees  of  freedom,  but  no  provisions  for  vibration, 
electronic  excitation,  or  nuclear  excitation.  At  moderate 
temperatures,  however,  both  quantum  mechanics  and  statistical 
thermodynamics  show  the  rigid  rotor  to  be  an  excellent  model  of 
the  diatomic  gas  molecule. 

Beginning  with  the  classical  Liouville  equation,  the  author 
develops  a  "Boltzmann-type  equation"  governing  the  motion  of  a 
rigid  rotor.  With  the  exception  of  four  additional  independent 
variables  due  to  angular  momentum  and  orientation,  the  Boltzmann- 
type  equation  is  identical  to  the  Boltzmann  equation.  The 
Boltzmann-type  equation  is  solved  via  the  Hilbert-Enskog-Chapman 

theory1'2  for  the  first  two  terms  of  the  expansion  of  the 
distribution  function.  The  governing  equations  of  motion  of  the 
gas  are  then  derived  by  substituting  those  first  two  terms  of  the 
distribution  function  expansion  into  the  conservation  equations 
of  mass,  linear  momentum,  total  angular  momentum,  and  energy. 

The  resulting  equations  of  motion  mimic  the  Navier-Stokes 
equations  except  with  regard  to  the  additional  angular  momentum 
conservation  equations,  the  molecular  constants,  the  heat  flux 
vector,  and  the  stress  tensor. 

The  development  of  the  Navier-Stokes  equations  from  kinetic 
theory  is  based  upon  a  smooth-sphere  molecular  model  possessing 
three  translational  degrees  of  freedom  only.  Such  a  model  is  an 
excellent  representation  of  monatomic  gas  molecules  at  moderate 
temperatures.  The  thermodynamic  properties,  i.e.  equilibrium 
properties,  of  monatomic  gases  differ  significantly  from  those  of 


M 


the  diatomic  gases.  Table  1  presents  this  difference  and  also 
the  excellent  comparison  between  experimental  results  and 
theoretical  calculations  based  upon  statistical  thermodynamics 
and  the  above  molecular  models.  One  purpose  of  this  study  is  to 
investigate  whether  weak  non-equilibrium  properties,  i.e.  flow 
properties,  also  strongly  depend  upon  the  gas  type  in  question. 

If  the  flow  properties  do  strongly  depend  on  the  gas  type,  the 
Navier-Stokes  equations  might  be  inaccurate  when  applied  to  some 
air  flows.  However,  we  shall  see  that  this  is  not  the  case. 

Another  purpose  of  the  current  work  is  to  introduce  greater 
rigor  into  the  development  of  the  Boltzmann-type  equation  than 

had  other  workers?'6'7  Kirkwood8'9  paved  the  way  in  this  regard; 
the  present  work  merely  specializes  his  development. 

Some  preliminary  definitions  and  conventions  are  in  order. 
The  present  work  deals  almost  exclusively  in  a  five-dimensional 
Euclidean  space.  A  "surface"  in  this  space  is  defined  to  be  a 
four-dimensional  geometry,  while  a  "volume"  is  a  five-dimensional 
geometry.  Both  vector  and  indicial  notations  are  used  in  the 
present  derivation.  The  specific  selection  of  one  notation  over 
another  is  based  upon  the  ease  of  interpretation.  By  convention, 
when  the  integrand  is  a  function  of  the  momentum,  the  integral 
limits  are  from  negative  infinity  to  positive  infinity  unless 
otherwise  stated. 

The  governing  equations  of  motion  are  next  derived. 

Important  results  and  notable  assumptions  are  then  summarized  and 
the  paper  closes  with  a  discussion  and  a  short  conclusion. 

Several  mathematical  developments  appear  in  the  appendices  at  the 
rear  of  the  paper. 


Table  1.  Molecular  Constants 

MONATOMIC  GASES 


exp 

calc 

exp 

calc 

exp 

calc 

V 

Y 

Pro°c 

Pr0°C 

S25°C 

S25°C 

He 

1.666 

1.667 

0.664 

0.667 

126.3 

125.3 

Ne 

1.64 

1.667 

0.675 

0.667 

146.5 

145.5 

Ar 

1.666 

1.667 

0.659 

0.667 

155.0 

154.0 

Kr 

1.68 

1.667 

0.656 

0.667 

164.2 

163.2 

Xe 

1.67 

1.667 

0.649 

0.667 

169.8 

168.9 

DIATOMIC  GASES 


exp 

calc 

exp 

calc  exp 

calc 

Y 

Y 

Pro°c 

Pr0°C  S25°C 

S25°C 

H2 

1.405 

1.400 

0.710 

-  130.7 

129.7 

HC1 

1.39 

1.400 

0.881* 

-  186.2 

186.0 

N2 

1.405 

1.400 

0.722 

-  191.5 

190.8 

CO 

1.402 

1.400 

0.745 

-  197.5 

196.8 

°2 

1.396 

1.400 

0.718 

-  205.1 

195.1 

C12 

1.35 

1.400 

0.771* 

-  223.0 

219.8 

* 

Experimental 

results 

at  25°C . 

Entropy  is  in  units  of  J7°K-mole.  Results  gathered  from  Refs. 
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SECTION  II 

DEVELOPMENT  OF  THE  GOVERNING  EQUATIONS  OF  MOTION 
1.  Derivation  of  the  Boltzmann-type  equation 

The  Liouville  equation2 


<£)’ r _ 

11“  , 

<0t 

Jt 

Lk  yi‘ 

i=i 

for  the  distribution 

f  unction 

Yf(NI<pN.3N,t 

space  (p^,qjj)  of  systems  of  N  molecules  is  the  fundamental 


equation  of  classical  statistical  mechanics.  The  Liouville 

v  ( N ) 

equation  cannot  generally  be  solved  for  Tf  .  Fortunately,  such 
detailed  information  as  that  contained  in 
v  ( N ) 

f  is  not  needed. 

The  present  analysis  shall  be  restricted  to  the  first-order 
distribution  function,  ^f(M  of  molecule  a  in  the  sub-phase  space 

ct 

(p  ,q  )  of  v-space.  Since, 
a  a 

.  If 

the  distribution  function  Yk  1  *  is  simply  the  integral  of  Yf(N) 

a 

over  the  sub  phase  space  (p^q^)  where  g  is  the  set  of  all 
molecules  exclusive  of  molecule  a: 

vr‘'U-  .x  JU‘“L  -  ,m.  u 


,(i\  CW) 


V..-'  .''.'"‘.''.■■'aV-.v.'-v. 


_ vV'/v  vlv.'Xvlv.  . 


provided  such  an  integral  is  finite.  This  requires  Yf  ' 


-►  0 


sufficiently  fast  as  |p|,  |q|  -»  «.  The  Liouville  equation  in 

terms  of  Tf  is  thus 
a 

.  jw  J'C  jj/C-o 

J3t  Jt  Jf*  Jf. 


The  Liouville  equation  is  next  specialized  to  the  case  of  a 
gas  composed  of  rigid  rotor  molecules  by  substituting  for  the 
Hamiltonian  of  a  rigid  rotor.  The  Hamiltonian  of  a  rigid  rotor 

in  an  inertial  reference  frame  is3: 


5 


where  x,  y,  and  z  are  the  Cartesian  coordinates  of  the  rotor's 
center  of  mass  and  0,  $  are  the  Euler  angles  of  the  rotor's 
orientation,  such  that 


(K0<ff 
01012  it. 

He  is  the  potential  due  to  effects  external  to  the  gas.  H1,  the 
intermolecular  potential,  is  a  function  solely  of  the  relative 
locations  and  orientations  of  the  molecules.  Substituting  the 
rigid  rotor  Hamiltonian  into  equation  four  yields: 


,’V  .T  ■■*  A  .'V'.V  JVJ.'  V,’.'  V’  U.*X*  ’-VTV'-.'V-V  ’ 


where 


1*1  X'  p(  3  mu.*,  B,  *  *n 

i=2  =%>  *}-  4vr*c\i»n  '  fxs  ***  j  Bg3  ** 

js3  "=^>  gr  *  •  p3  r  mw  •  8 j  =  m 

i=  ^  ^  ©  -  <i.V*c4l6v,  ;  =■  _L©  j  Bif  *  X 

j  =  5  =£>  <p  -  <l,mW  ;  j»5  r(l  siw*e)g>  j  B^-  :  T  S»*»*© 


The  last  terra  of  equation  six  represents  a  torque  causing  the 
rotor  to  precess  from  a  general  rotation  to  a  rotation  where  9  is 
fixed  at  ir/2.  However,  since  the  momental  ellipsoid  of  a  rotor 
is  an  infinitely  long  circular  cylinder,  the  free  rotation  of  a 
rigid  rotor  is  stable.  Therefore,  this  precessional  term  exists 
only  in  the  presence  of  external  torques.  Such  torques  arise 

from  both  He  and  H1.  The  present  analysis  shall  be  restricted  to 
nonpolar  molecules  undergoing  impulsive  collisions  so  that  the 
precessional  term  vanishes.  Since  nonpolar  molecules  are 
unaffected  by  external  field  torques,  however,  forces: 


rife. 


*  O. 


The  above  two  assumptions  lead  to  the  equation: 


\*rw 

Z2*.  ,  Pi*. 


Since 


7 


8 


9 


v.;< 


y'\'v 


equation  eight  may  be  rewritten  with  the  aid  of  equation  two  as: 


The  right-hand  term  of  equation  10  is  the  collision  integral. 

The  collisional  integral  must  be  reduced  to  be  a  function  of 

(  1  )  8  9 

f  •  Appendix  A  follows  Kirkwood's  work  '  for  the  reduction. 

By  assuming  the  gas  to  be  dilute  and  in  molecular  chaos,  we  find 
equation  10  to  be  equal  to: 


+  Ik. 


where : 


and  where  x  is  a  characteristic  time  long  in  comparison  to  the 

t 

duration  of  a  collision  but  short  in  comparison  to  a  Poincare 
cycle8 . 

Appendix  B  develops  the  actual  equations  of  motion  governing 
binary  collisions  between  rigid  rotors.  For  any  given  direct 
collision  between  impulsively  colliding  rigid  rotors,  we  see  that 
there  exists  an  inverse  collision.  The  existence  of  inverse 
collisions  is  of  crucial  importance  to  the  present  work. 

Boltzmann's  H-Theorem1  may  be  shown  to  hold  for  rigid  rotors  and 
the  relation 


dpadpb-  dp^dp^ 


holds . 


The  relative  spatial  coordinate,  qR,  deserves  further 
consideration.  The  coordinate  is  a  five-space  volume.  Only  the 


-v-\ 


portion  of  five  space  in  which  a  type  b  molecule  is  able  to 
collide  with  a  type  a  molecule  within  time  t  contributes  to  the 
collision  integral.  This  portion  of  five  space  may  be  considered 

as  a  five-space  cylinder  with  the  generator  |\?R|.  The  vector, 

\L,  is  the  relative  velocity  of  the  points  of  collision.  From 

K 

the  development  in  Appendix  C,  it  is  straightforward  to  show 

cbj^=  iVjjl f  sihCsi^T  dr  icd X 

is  the  volume  of  the  five-space  cylinder  in  question.  Equation 
11  is  thus: 

it  <U  .  15 

Special  attention  must  be  paid  to  the  regularity  of  the  base 
of  the  five-space  cylinder: 


jL  T  fir  ^  - 

Jt  J*}u.  ‘Vi* 


This  surface  may  not  be  regular  for  a  true  rotor.  In  terms  of 
the  variables  of  Appendix  B,  the  surface  is  piecewise  continuous 
when 

cosy  ill!!  17 

Therefore,  the  regularity  of  the  collision  surface  is  in  question 
until  a  transformation  is  made  to  from  the  variables  of  Appendix 
B  to  the  variables  of  Appendix  C.  If  the  collision  surface  is 
not  regular  along  the  curve  in  question,  the  collision  integral 
will  have  to  be  separately  evaluated  over  each  regular  surface. 


Equation  15  is  identical  in  form  to  the  Boltzmann  equation 
but  has  an  additional  four  independent  variables  corresponding  to 
the  orientation  and  angular  momentum  of  the  rotor.  Due  to  this 
similarity,  equation  15  shall  be  called  a  Boltzmann- type 
equation. 

2.  Derivation  of  the  Conservation  Equations  of  a  Rigid  Rotor 

The  expansion  procedure  of  Enskog  may  be  extended  to  solve 
the  Boltzmann- type  equation.  We  begin  be  developing  the 
conservation  equations  for  a  gas  composed  of  rigid  rotor 
molecules . 

Let  ¥ ( p  , q  , t )  be  any  scalar,  vector,  or  tensor  function 

OL  GL 

of  the  molecular  momentum,  position,  and  time.  The  change  of  ¥ 
per  unit  time  per  unit  volume  at  the  point  qQ  due  to  collisions 
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[>(^  ip 

J  rUt  L*b 


where  n  is  the  number  density  such  that 


j-HV  N 


By  equation  15,  we  have: 


where  the  right-hand  side  of  equation  20  is  the  change  of  ¥  per 
unit  time  at  the  point  q  due  to  molecular  flow.  By  simple 
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manipulations  (see  Appendix  E) ,  equation  20  may  be  transformed 
into  Enskog's  general  equation  of  change  of  the  property 

¥  ( pa ,  qa , t )  : 
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Note  that  the  number  average  of  the  molecular  properties,  ®ja' 
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Unless  noted  otherwise,  B.  will  now  denote  B.  . 

]a  ja 

For  properties  which  are  conserved  during  collisions 

at  =  o. 
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The  conserved  properties  are  the  number  of  molecules  (or  mass), 
linear  momentum,  total  angular  momentum,  and  energy.  Thus: 
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Upon  decomposing  the  momentum  into  a  mean  component: 


and  into  a  peculiar  component 
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and  substituting  for  each  conserved  property  into  equation  22 
yields  the  four  conservation  equations: 

i)  Continuity  Equation: 


*  -hit  *  i 
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ii)  Conservation  of  Linear  Momentum: 
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iii)  Conservation  of  Total  Angular  Momentum: 
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iv)  Conservation  of  Energy: 
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27 


28 


=  O  29 


where  the  subscripts  have  been  deleted  for  notational  ease. 
Equations  26,  27,  28,  and  29  are  fully  derived  in  Appendix  E. 
Equations  27  and  28  show  that  conserving  linear  momentum  and 
total  angular  momentum  at  a  point  is  equivalent  to  separately 
conserving  linear  momentum  and  internal  angular  momentum.  We  may 
thus  combine  equations  27  and  28  to  yield: 


v)  Conservation  of  Total  Momentum  at  a  Point: 
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The  five-space  stress  tensor  and  the  heat  flux  five  vector 
are  by  definition: 
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respectively.  Both  . 

and  Qj 

are 

determined  from  the 

distribution  function,  f 


3.  Enskog's  Expansion  Method  of  Solution 

Under  Enskog's  expansion  method  of  solution,  the 
distribution  function  is  assumed  expandable  into  the  convergent 
infinite  series: 
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The  present  analysis  derives  the  first  two  terms  of  the  series 
f1  and  f°. 

Let  the  Boltzmann-type  equation  be  defined  as1: 


where 


Enskog  supposed  that  when  the  Boltzmann-type  operator,  £, 

—  f  V\ 

operates  on  the  above  infinite  series  expansion  of  f,  the  r— 
term  of  the  resulting  series  depends  only  on  the  first  r  terms 


s  $)  -  srr*  ?■  >v ■* ...)  -  rrn  *  m'.n  *  s\ > ... 


The  functions  f  are  also  assumed  to  satisfy  the  separate 
conditions : 


rdVo 

... 


Enskog  subdivided  £  such  that 


S'- O'*  r-o  w  ^1. 
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Equation  40  may  be  recast  as: 

The  general  solution  to  equation  41  is  simply  the  sum  of  the 
particular  solution  and  the  homogeneous  solution.  Let  the 
homogeneous  solution  be 

r-TT 


where  r  is  the  homogeneous  solution  to  equation  41.  Then: 
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and  therefore  9  must  be  a  linear  combination  of  the  conserved 
properties  of  a  collision,  i.e. 

Xr_  M  V  !*•  x 


lr  *• r  V  !*•  x  , 
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where  ocl'r  and  a 4,r  are  scalars,  a2,r  is  a  three-vector,  and  a3,r 
is  a  five-vector.  The  above  constants  are  chosen  such  that 


Irrv 


for  i=l,2,3,4;  r>l 


which  ensures  that  the  average  values  of  n,  T,  and  pQ  are  given 


jrfij 


for  i=l,2,3,4. 
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The  preceeding  brief  overview  of  Enskog's  method  of  solution 
brings  us  to  the  important  subject  of  the  conditions  of 
solvability  of  equation  40.  If  equation  40  is  satisfied,  then 


for  i=l,2,3,4, 


By  the  developments  of  Ref.  1,  J  is  divisible  into  pairs  of 
terms  such  as  J(fpf^)  +  J(fqfp>  plus  a  term  J(fpfp)  if  p  is  even. 
It  is  straightforward  to  show  that  the  integral  of  the  product  of 
•F1  and  these  terms  vanishes  because 


during  a  collision, 
of  equation  40  is: 


¥  +  *b  «  1 "  + 


Thus,  a  necessary  condition  for  the  solution 
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for  i=l,2,3,4. 


Hilbert1  proved  that  the  above  condition  is  also  a  sufficient 

condition  for  the  solution  of  equation  40  by  showing  I(4r)  may  be 
expressed  as 


.  j  -n.tp.p0  i*Pt 


where  ( p , p^ )  is  a  symmetric  function  of  p,  p^.  Such  a  proof  is 

not  yet  available  for  the  present  case.  However,  it  is  fully 
expected  that 
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where  dqR  is  given  by  equation  30  is  indeed  expressible  in  the 
form  of  equation  49.  By  following  arguments  identical  to  those 
of  Chapman  &  Cowling1,  it  is  straightforward  to  show  that 

\\  ■  -a.tfilVp'i .  si 

It  remains  to  show  that 

where  il(p,p^)  is  symmetric  with  respect  to  p,  pfe.  This  is 
strongly  felt  to  be  true.  It  will  be  shown  that  f°f°  i3 
symmetric  with  respect  to  p,  p^.  Since  inverse  collisions  exist. 


the  problem  reduces  to  whether  or  not  dq^  is  symmetric  with 

■4  — ►  — ♦  -4 

respect  to  p,  pfc.  If  dqR  is  not  symmetric  in  p,  pfa,  then  the 

principle  of  objectivity  would  be  violated,  i.e.  the  final  answer 
would  depend  upon  one's  choice  of  the  origin  and  coordinate 
system.  Physically,  this  does  not  happen  and  one  would  therefore 

conclude  that  dq^  is  indeed  symmetric  in  p,  p^.  A  mathematical 

proof  is  required,  however,  to  rid  the  present  work  of  any 
ambiguities  in  this  area.  The  present  work  shall  proceed  on  the 
reasonable  assumption  that  equation  48  is  a  sufficient  condition 
for  the  solution  of  equation  40. 


4.  Derivation  of  f° 


The  first  term  f°  of  the  series  expansion  of  f  is  easily 
found. By  Enskog's  method  of  division, <D* =0 ,  and  the  Boltzmann- 
type  equation  is: 

O.  j[[(j  (VX  -  TOW,)  Wi  i.  Jp  it  U  . 

Thus,  in  general1 
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or 
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Equation  55  can  hold  in  general  only  if  ln|f  |  is  a  linear 
combination  of  the  properties  conserved  during  collision.  These 

properties  are  ¥1  for  i  =  1,  2,  3,  4, 


Thus : 


•  r  ^  ^  l.t>  m*  <4  ft  1  . 
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By  the  development  of  Appendix  F: 
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where  B5  is  Isin  0,  not  1/2.  The  equilibrium  distribution 
function  is  built  on  the  following  definitions: 

f  =  \n 
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5.  Derivation  of  f 

The  evaluation  of  the  second  term,  f1,  is  more  difficult. 

By  Enskog's  method  of  division 

_$'  =  &  J‘  *  O 

where  £>'  operates  on  f°  only.  Since  f°  is  a  function  only  of  the 
peculiar  momentum,  we  may  recast  <D'  in  terms  of  P  as: 

w.  ii.J  j*-  .  h. jl* .  h  &  j*-  .  62 

The  derivation  of  equation  62  is  not  presented  as  it  is  identical 

to  the  derivation  presented  by  Chapman  &  Cowling1.  The 

equilibrium  distribution,  f°,  is  already  known  and  the  first 
approximations 

E.Pol  D«T 
Dt  ’  Dt 

may  be  found  from  the  conservation  equations  and  Enskog's 
subdivision  of  <D  .  The  conservation  equations,  equations  26,  30, 
and  29,  may  be  recast  as: 
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respectively.  Enskog  divided  these  time  derivatives  as  follows 
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where  A  may  be  either  n,  pQ ,  or  T.  The  division  of  the  right- 


hand  side  is: 
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By  defining  the  first  approximation  to  the  substantive  derivative 
as : 

^2  £  *  P«a  jsL 

Dt  Jt  6i 


equations  62  and  63  become: 
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When  r  is  zero,  only  the  even  terms  of  the  integrands  of 
equations  74  and  75  survive  the  integration.  Thus, 
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77 
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=»nd  by  the  definitions  and  results  of  Appendix  F: 
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By  defining  the  thermodynamic  pressure  as 


.•itn  •  mff«  -  -  2  ^ 


The  thermodynamic  pressure  may  not  equal  the  hydrostatic 
pressure.  However,  the  hydrostatic  pressure  cannot  be  deduced 
because  Maxwell's  relations  are  strictly  valid  for  a  three- 
dimensional  space  only.  Thus,  Maxwell's  relation: 


p  s  _  ( — )  =  kT  ( ) 
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cannot  be  used  in  the  present  study. 

Substituting  the  first  approximations  to  the  heat  flux 
vector  and  the  stress  tensor  into  equations  76  and  77  yields: 
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Substituting  equation  82  into  equation  62  yields: 
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Appendix  H  reduces  equation  84  and  finds  that  £1=0  is  equivalent 


where  ,  the  reduced  momentum,  is  defined  as: 
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Appendix  H  derives  the  form  of  $  as: 
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where  J((?)  and  K((?)  are  scalar  functions  of  n,  and  T. 
Substituting 

V-iV 

for  f  in  the  definitions  of  the  five-space  stress  tensor  and  the 
heat  flux  five  vector  determines  the  second  order  approximations 

of  and  Qy  Appendix  I  details  the  derivations  of  and  . 

The  resulting  expressions  are: 


=  -2 


where 


m.  = 


tr 

2 ? 


88 


89 


90 


and 


91 


92 


22 


f 

%*  \* 


$3 


£ 


>> 

>.v 
■  * 


w-V. 

W 

vV 

atz-ji?*? 


«/ 

>>> 

■'■'i 

v  v 

v.y; 

V:W 

7',.’ 

f.O 

f.V.'v 

f’.S.'-' 

*.V.\ 


v> 

l"  v  *1 

V>v 

‘."-NV 


»*1 

~:”V_ 


.  v.v. 


N^v; 

V  •/ 
\  \  ■ 

- 

v'? 


V*  Aw".  V*  * 


*  '  V  •,  • 


where 
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Appendix  J  evaluates  y  and  A  through  the  use  of  Sonine 
polynomials.  Assuming  that  the  first  term  of  the  Sonine 
polynomial  is  of  sufficient  accuracy.  Appendix  J  shows: 

A’  ££  fe* I 


where  the  terms  a^  and  b^  can  be  evaluated  only  when  the 

intermolecular  potential  is  known.  Since  the  intermolecular 
potential  of  the  present  problem  is  simply  a  pair  of  Heaviside 
step  functions,  a^  and  b^  could  be  evaluated.  Such  an 

evaluation  is  likely  to  be  very  difficult.  Given  the  assumption 
of  impulsive  collisions,  the  evaluation  is  also  prone  to  be  in 
significant  error  with  respect  to  the  experimentally  measured 
values.  Therefore,  neither  a^  nor  b^  have  been  evaluated. 
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SECTION  III 


SUMMARY  OF  DERIVED  EQUATIONS,  CONSTANTS,  AND  ASSUMPTIONS 


The  equations  of  motion  of  a  diatomic  gas  whose  molecules 
may  be  modelled  as  rigid  rotors  are: 


1)  Conservation  of  the  Number  of  Molecules: 


ia  *  fell  ~ 


2)  Conservation  of  Total  Momentum  at  a  Point: 


W  i.i s 


where : 


3)  Conservation  of  Energy: 
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where 
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4)  Equation  of  State: 

pr*HT 

where  p  is,  by  definition,  the  thermodynamic  pressure. 
Additional  quantities  of  interest  are: 

i)  Specific  Heats: 

=  ;  Cp=3.S„k 

ii)  Ratio  of  Specific  Heats: 

y  .  1.100 

iii)  Prandtl  Number: 


iv)  Equilibrium  Distribution  Function: 

.  *  TT  few  C*  f? / 2 

irl 

Most  assumptions  made  in  the  present  analysis  are  identical 
to  those  of  the  derivation  of  the  conventional  Navier-Stokes 
equations.  Those  assumptions  unique  to  the  present  analysis  are 
listed  below: 


1)  The  molecules  nonpolar.  This  amounts  to  setting  F®  =  F^  =  0. 

2)  The  intermolecular  collisions  are  impulsive. 

3)  That 

I  r  o 


is  a  sufficient  condition  for  the  Enskog  expansion  method  of 
solution  of  equation  15.  This  appears  to  be  a  reasonable 
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assumption  based  on  physical  arguments,  the  form  of  f°,  and  the 
existence  of  inverse  collisions.  Strictly  speaking,  however,  the 
sufficiency  condition  is  an  assumption  and  a  rigorous  proof  is 
needed. 
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SECTION  IV 


DISCUSSION 

We  have  seen  that  a  nonpolar  diatomic  gas  whose  molecules 
collide  impulsively  obeys  a  Boltzmann-type  equation,  namely 
equation  15.  The  differential  portion  of  the  equation  is  easily 
developed  by  substituting  the  well  known  Hamiltonian  of  a  rigid 
rotor  into  the  Liouville  equation.  The  collision  integral  is 

8  9 

developed  along  the  lines  of  Kirkwood's  work.'  Appendix  B 
develops  the  actual  governing  equations  of  motion  of  a  binary 
collision  between  rigid  rotors.  Appendix  B  displays  the 
fundamental  point  that  inverse  collisions  exist  as  long  as  the 
collision  mechanics  are  consistent  with  the  Hamiltonian.  The 

rough  sphere  molecular  model6  suffers  in  this  regard,  as  the 
Hamiltonian  does  not  account  for  the  strain  energy  of  the  rough 
sphere  during  collision. 

The  derivation  of  the  differential  portion  of  the  Boltzmann- 
type  equation  is  not  the  straightforward  extension  from  three 

degrees  of  freedom  assumed  by  others?'7'  In  fact,  an  additional 
term  describing  the  precession  of  the  rotor  is  generally  present. 
Only  when  the  molecules  are  nonpolar  and  collide  impulsively  does 
this  precessional  term  vanish.  Furthermore,  it  is  unlikely  that 
a  general  Boltzmann-type  equation  can  be  derived  for  a  rigid,  six 
degree  of  freedom  molecule.  Consider  the  classical  Hamiltonian 

for  such  a  molecule3: 

H  =  2^Cp x  +py  +  pz3  +  ~~"T 2Q[(p4»  “  P*cose>cosiK  -  p^inesimj,] 2 

21^sin  6 

+  ;T  1  2Q[(p$  "  P*cos0>sin*  +  p0sin0cosip] 2  +  PJ 
Z ln s in  0  C 


where  I^,  Ig,  and  1^,  are  the  principal  moments  of  inertia;  0,  0, 
and  0  are  the  Euler  angles,  and 

O<0<.ir 

Oi.0i.2ir 

0i.iJ)i.2ir. 

Substituting  the  above  Hamiltonian  into  the  Liouville  equation 
produces  a  complicated  equation  in  f.  There  are  many  terms 
describing  the  precession  of  a  general  rotation  to  a  rotation 
about  an  axis  of  minimum  or  maximum  moment  of  inertia.  These 
precessional  terms  do  not  disappear  when  the  molecule  is  assumed 
to  be  nonpolar  and  to  collide  impulsively.  The  rigid  rotor 
molecule  thus  appears  to  be  the  most  complex  molecule  that  obeys 
a  Boltzmann -type  equation  with  a  minimum  of  assumptions. 

Molecular  models  of  six  degrees  of  freedom  and  more,  such  as  the 
non-rigid  rotor,  the  symmetric  top,  and  the  asymmetric  top,  may 
require  an  approach  based  upon  a  pair  or  triple  distribution 
function. 

The  assumptions  of  the  present  study  mirror  those  in  the 
development  of  the  Navier-Stokes  equations,  except  in  four 
instances.  One  of  the  assumptions  different  from  those  of  the 
Navier-Stokes  equations  is  the  basis  of  the  present  study,  namely 
that  the  gas  molecule  can  be  modelled  as  a  rigid  rotor.  The  two 
additional  assumptions  of  impulsive  collisions  and  nonpolar 
molecules  eliminate  the  precessional  term  from  the  Boltzmann-type 
equation.  The  fourth  assumption  that 

is  a  sufficient  condition  for  the  solution  of 


is  critical  to  the  development  of  the  stress  tensor  and  the  heat 

flux  vector.  If  that  assumption  is  false,  then  I(*r)  is 

asymmetric  with  respect  to  pa  and  p^.  The  collisional  equations 

of  motion  would  then  change  upon  an  interchange  of  the  roles  of 
molecules  a  and  b.  The  physics  of  the  motion,  however,  does  not 
depend  upon  how  we  select  molecules  a  and  b.  Therefore,  this 
sufficiency  assumption  is  probably  correct.  A  rigorous  proof 

similar  to  Hilbert's1  treatment  in  three  degrees  of  freedom  is 
required,  however,  to  close  the  discussion  on  this  issue. 

The  stress  tensor  is  expanded  in  Fig.  1.  Those  terms 

containing  angular  quantities,  0,  0,  10,  and  10,  are  periodic 
because  the  Euler  angles  are  periodic.  Also,  there  are  seven 
coefficients  of  “viscosity".  These  coefficients  of  viscosity  and 
their  orders  of  magnitude  relative  to  the  usual  shear  viscosity, 
m,  are  displayed  in  Table  2. 

Table  2.  Relative  Magnitudes  of  the  Viscosity  Coefficients 
Coefficient  of  Viscosity  Relative  Order  of  Magnitude 


ress  Tensor 


The  molecular  radius  of  gyration  is  r  .  The  relative  orders  of 
magnitude  of  the  seven  coefficients  of  viscosity  are  based  upon  a 
typical  value  of  r  of  10”^  meters.  The  first  three  columns  of 

g 

the  stress  tensor  generate  changes  in  linear  momentum.  We  see 
the  periodic  terms  probably  contribute  negligible  changes  in 
linear  momentum  because  of  their  very  small  coefficient  of 
viscosity.  The  periodic  terms  will  perhaps  be  significant  only 
in  shock  layers  where  the  gradients  are  large. 

One  should  also  note  that  a  gas  composed  of  rigid  rotor 
molecules  may  have  a  bulk  viscosity.  The  bulk  viscosity  will  be 
due  to  the  rotational  degrees  of  freedom  responding  to  energy 
gradients  at  a  rate  different  from  the  response  rate  of  the 
translational  degrees  of  freedom. 

The  derivation  of  the  Prandtl  number  is  unfinished.  The 
ratio  of  is  unity  for  the  rigid,  smooth-sphere  molecular 

model.  If  this  is  found  to  hold  for  rigid  rotors,  the  Prandtl 
number  will  be  .8  --  a  mediocre  prediction.  See  Table  1. 

The  present  work  forms  the  foundation  upon  which  a  study 
into  the  properties  of  a  gas  composed  of  non-rigid  rotors  can  be 
built.  It  would  be  most  interesting  to  establish  the  Liouville 
equation  in  terms  of  a  pair-distribution  function  and  then 
substitute  for  the  Hamiltonian  of  each  atom  of  the  rotor.  The 
Hamiltonian  must  include  the  quantized  intermolecular  potential 
that  holds  together  the  rotor.  If  this  approach  is  possible  not 
only  vibration,  but  disassociation,  could  be  included  into  the 
governing  equations  of  motion.  The  resulting  equations  should 
model  the  important  real  gas  effects  of  hypersonic  flows. 


SECTION  V 


CONCLUSION 

One  may  develop  the  equations  of  motion  of  a  diatomic  gas  by 
substituting  the  Hamiltonian  of  a  rigid  rotor  into  the  Liouville 
equation.  The  equations  mirror  the  Navier-Stokes  equations 
except  with  regard  to  an  additional  two  equations  conserving 
angular  momentum,  the  stress  tensor,  and  certain  molecular 
constants.  The  assumptions  of  the  present  derivation  are 
identical  to  those  in  the  development  of  the  Navier-Stokes 
equations  except  that  the  molecules  are  nonpolar  and  collide 
impulsively.  Intriguing  because  it  resides  in  a  five-dimensional 
space,  the  development  is  a  stepping-stone  to  a  study  which  will 
include  the  vibration  and  disassociation  of  a  diatomic  gas. 


By  Gibbs'  postulate,  the  ensemble  average  of  any  molecular 
function,  is  the  roacroscopically  observed  value: 

Xws  ’  8  ^  A 


0 

Kirkwood  argues  further  that  the  measurement  of  any  macroscopic 
property  is  an  inherent  time  averaging  operation.  Therefore, 

A-2 

where  the  time  interval,  t,  is  long  in  comparison  to  collision 

durations  but  short  in  comparison  to  a  Poincar6  cycle8.  Kirkwood 
also  shows  that  the  measured  value  of  any  molecular  function, 

^a'V  is 


A. 3 


where 


A. 4 


Recasting  equation  1 1  of  the  text  in  terms  of  the  time-averaged 
distribution  function  leads  to: 


'sr.  r:  r.T,  v.’ttct: ?■  ~ /  "? 
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Through  the  use  of  the  Dirac  delta  function,  the  collision 
integral  may  be  extended  over  all  -y- space: 


A. 


*  iKi .  ♦  A  J  ^  s 

where 

*1»1  t;is' 

The  subscript  o  denotes  a  value  at  some  initial  instant  of  time, 
t.  By  the  principle  of  conservation  of  extension  in  phase  space 


A. 


The  gas  shall  be  considered  sufficiently  dilute  30  that 
binary  collisions  are  the  sole  mechanism  of  the  collisional 
transfer  of  molecular  properties.  Since  the  collisions  are 
impulsive,  the  interraolecular  potential  is  a  set  of  Heaviside 
step  functions.  The  potential  due  to  any  molecule  a  is  zero 
outside  a  certain  volume  of  influence,  V.  Therefore,  in  the 
region  of  configuration  space  exclusive  of 

H 

Lvl 


the  gas  is  in  "molecular  chaos".  Kirkwood  shows  that  for  a 
dilute  gas  in  molecular  chaos: 


,o'C  d 


[‘i,-  - 


8 

$8 

& 


A. 


iffwv* 


By  way  of  equation  two  of  the  text,  this  leads  to: 
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Since  the  intermolecular  collisions  are  impulsive,  the  time 
interval  t  may  be  made  arbitrarily  short.  Thus: 

H  * 1  v  V7.  V 

'tJ>o  J 


A. 10 


Assuming  the  above  to  be  true  for 


equation  A. 9  may  be  integrated  at  once: 

?■  fef 

By  Newton's  second  law  of  motion,  KirJcwood  recast  the  collision 
integral  as: 


jfjf)  t  ! 


A. 12 


A. 13 


The  gas  shall  be  considered  to  be  composed  of  identical 
indistinguishable  molecules.  Equations  A. 12  and  A. 13  may  be 
recast  to  apply  to  molecule  phase  space,  p-space.  Noting: 
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A. 14 
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A. 15 


leads  to: 


A. 17 
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A. 18 
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Without  any  loss  of  generality,  q^Q  may  be  expressed  relative  to 

qaQ.  Denoting  this  relative  spatial  coordinate  qRo  and  dropping 
the  superscripts  from  the  p-space  probability  density  leads  to: 
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This  may  be  integrated  at  once: 
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since,  by  Newton's  third  law  of  motion: 


By  invoicing  the  principle  of  conservation  of  extension  in  phase 
space  and  by  performing  a  second  time-averaging  operation  over  x 
gives 

•  i* £  ■  ^  A- 21 


It  can  be  shown  that 


A. 22 


As  mentioned  previously,  that  the  collisions  are  impulsive  allows 


Thus : 
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Without  loss  of  generality,  we  may  locate  the  origin  of  our 
coordinate  system  at  the  c.g.  of  molecule  a  at  the  instant  of 

impact.  Further,  we  may  also  rotate  the  axis  system  to  align  x_ 

Q 

with  the  axis  of  molecule  a.  Thus: 


The  collision  geometry  of  molecule  b  with  respect  to 
molecule  a  i3  uniquely  described  by  the  quantities: 

c 
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ft, 

-A 

where  Tfc.  tx--  Irjcoiy  . 
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The  equations  of  motion  are: 
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Assume : 
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The  above  assumptions  in  conjunction  with  equation  B.5 
yields 


B. 

-*  1  -*  -k 

^  vb*"  X 

B. 

X^a  "X^a*  ro.X-J" 

B. 

IwJ  -  lAj"  rb*X 

B. 

There  are  four  equations,  B.6  through  B.9  in  five  unknowns,  7, 
V  ,  ,  u.  ,  w.  .  The  fifth  equation  of  motion  depends  upon  the 
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rod  material  and  surface.  The  rods  shall  be  assumed  to  be 
perfectly  elastic  and  perfectly  smooth.  Thus,  the  components  of 
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the  relative  velocity  tangent  to  the  "collision  plane"  are 
unchanged  by  the  collision  while  the  normal  component  is 
completely  reversed. 

The  plane  of  collision  is  that  plane  containing  the  three 


vectors  x  ,  r  and  r.  .  Let  the  normal  unit  vector  of  this  plane 


be  denoted  as  i 


N 


The  two  tangential  unit  vectors  will  be 


denoted  tT1  and  and  where  ^T1  is  oriented  such  that 

A  A 

ifpT •  i  =  1.  Thus,  the  final  equations  of  motion  are: 
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B.  11 
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B.  12 

and  the  relative  velocities  before  and  after  collision  are: 
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One  can  now  solve  for  J,  and  then  V  ,  V,  ,  w  and  u>,  .  By 
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equations  B.6  through  B.9: 
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By  equation  B.13; 
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Substituting  equation  B.15  into  equations  B.ll  and  B.12  shows: 
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Equations  B.16  and  B.17  can  hold  in  general  only  if  I1‘T\  andTli.T2 
which  also  states  "51?#.  and  Tl7b  .  Therefore,  tJIln  .  We  have  found 

the  direction  of  ?;  now  we  must  find  its  magnitude.  Substituting 
equation  B.15  into  equation  B.10  yields: 

j^R.'  ^  J^4  •  <-w  -  -V^-uN 
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Since  r  •  J  =  r.  ■  J  =  0: 
a  b 
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Since  I'Ln*±1T|  ,  define:  which  produces: 
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Thus  : 
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Substitution  of  equation  B.19  into  equations  B.6  through  B.9 
yields : 
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The  above  could  just  as  easily  be  cast  in  forms  where  the 
primed  variables  are  the  initial  quantities.  From  equation  B.10 
and  the  definition  that  Ve-iN  be  positive  in  3ign,  the 
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direction  of  i„  is  reversed  for  this  case.  Thus: 
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Therefore,  given  a  direct  encounter  where  the  initial  velocities 


are  ^ wa,  and  the  unit  normal  of  the  collision  plane  is 

i^,  there  exists  an  inverse  encounter  with  initial  velocities 

equal  to  the  final  velocity  of  the  direct  encounter  and  with  a 
unit  normal  vector  -iN  of  the  collision  plane.  Other  than 

reversal  of  the  collision  plane  unit  normal  vector,  the  geometry 

*4  -4  -4  -4 

of  the  two  collisions  is  identical  in  terms  of  r,  r  ,  rK,  and  z . 
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By  the  method  of  successive  projection: 
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Consider  the  integral 


J  t 

where  <?  is  the  five  vector 

C  =  aVVVWVRVS:*  1 

of  magnitude  C.  Let  0(C)  =  F(C)  U2.  By  symmetry: 

ju’fCc^c-jvVfOJc^^fHc^c  .  jaVCcW  *jffCO<l£  . 

-«•  -i  ‘t 

Then, 

ju2F(c)<lc  •  £ OcuNvW* R^S^fCcUc.  •  -Lj^FCcWc 
-*  *•  ] 

and 

j fcdctic  --  j  Ij'c'FCcUc. 

-A  -< 

because  the  integrals  of  the  non-diagonal  terms  vanish  as  they 
are  odd  functions  of  U,  V,  W,  R,  or  S.  Thus 

\  Kc\tl  4c  •  O  .  ] 

— c 

For  any  vector 

CfccxvcM^  *  K-fmiai  i 

.  ±Afe2r- cc^c.  ’] 

Next,  let  0(C)  =  S4F(C).  Where  S  =  Ccos3. 


Then: 
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Integrating  with  respect  to 
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and  next  with  respect  to  oc: 
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Thus : 
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The  integral 
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j  KaS*  F(C^C  *  FCOC^C 
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where  B  J*  BCc.^  may  next  be  evaluated. 
For  a  typical  diagonal  element: 
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while  for  a  typical  non-diagonal  element 
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Thus : 
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From  the  Liouville  equation 


*  JL .  a  &  , ;  li  .  ( 'iL\ 
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E. 


where  f  B.(q.),  is  generalized  f orce/torque ,  and 
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8, 

y-Z  ■» 
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-T 
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Let  ¥(p,q,t)  be  any  function  of  the  molecular  momentum,  position 
and  time.  The  function  may  be  a  scalar,  vector,  or  a  tensor 
function.  The  rate  of  change  of  ¥  per  unit  "volume"  due  to 
collisions  is: 


But : 


3  v\ 


E. 


and ,  expand i ng : 


since  p^pC^, 


E.  5 


Giving : 


*  A1 
E*  V 


E.  6 


If  ¥  happens  to  be  a  conserved  quantity,  then  A*F  =  0 .  Such 
quantities  are  mass  (for  non-reacting  gases  and  neglecting 
nuclear  effects),  linear  momentum,  total  angular  momentum,  and 
energy  (for  non-nuclear  effects).  These  conservation  equations 
f  ollow. 

i)  Conservation  of  Mass: 

Under  the  above  stipulations,  this  reduces  to  the  conservation  of 
the  number  of  molecules.  Thus,  ¥  =  1. 

Equation  E.6  becomes: 


4-  +  r  i^l i  =  O  •  E.7 
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Bs  is  now  1/2  for  the  above  equation  as  this  is  the  number 


average  value  of  13  .  Next,  decompose  the  momentum,  p,  into  a 

5 

peculiar  component,  and  a  mean  component  pQ .  Also, 


?  =  ?(p,q,t).  Therefore: 


?  =  h* P 


Equation  E.4  is  thus: 


A  *  I 


ii)  Conservation  of  Linear  Momentum: 


For  i  =  1,2,3,  =  pi  ,  whereas  ^  =  0  for  i  =  4,5 


Equation  E.6  yields: 
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By  continuity,  however. 
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iii)  Conservation  of  Total  Angular  Momentum: 
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Denote  ¥.  as  1.  for  k  =  1  -+  5.  Equation  E.6  yields: 
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E.  15 
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which  the  conservation  of  linear  momentum  crossed  with  the 
position  vector.  Thus,  for  i  =  1,2,3,  the  conservation  of  total 
angular  momentum  equation  is  redundant. 

For  i  =  4,5,  1^  =  and  by  the  development  of  the 

conservation  of  linear  momentum  equation: 


,  ..  -fry  Ifct  1 1  f  *  jljj  a  Q 


where  one  must  remember  i  =  4,5, 


iv)  Conservation  of  Energy: 


■t-  pj/«i 


Equation  E.6  yields: 
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Also : 


3  2f.ip.ji  *  ♦  ft-fj 


where 
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E.  32 


By  the  conservation  of  linear  momentum  and  total  angular  momentum 
equations: 
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Consider  equation  E.33  first. 
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By  the  conservation  of  mass,  the  above  reduces  to: 
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But  flu.  I  Ml  =  5tr  .  Thus 
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By  the  assumption  of  molecular  chaos  =  0.  Expanding  the  last 
term  yields 


Thus  : 
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Again,  by  the  assumption  of  molecular  chaos 


0.  In 


combination  with  the  conservation  of  mass  equations,  this  yields: 
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Next,  consider  equation  E.34, 


Bi  Jiilwi  ^Vis:)'  Si^iLi 


By  the  assumption  of  molecular  chaos,  this  reduces  to 
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Expanding: 
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By  molecular  chaos,  we  have: 
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Adding  equations  E.44  and  E.48  yields: 
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which  by  the  conservation  of  linear  momentum  and  total  angular 
momentum  yields: 
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By  Enskog's  expansion  method 


J3‘ 


=  o 


which  leads  to 


where 


F.l 


F.  2 


The  only  difference  between  the  above  collision  integral  and  the 
collision  integral  for  spherical  molecules  is  the  dimensional 
rank  of  the  vectors  ^  and  ^  .  Boltzmann's  H-theorem  is  rank 
independent  and  therefore  holds  for  equation  F.2.  Thus, 


or 


F.  3 


F.  4 


Since  lnf°  is  both  a  scalar  and  a  conserved  property  of  the  gas, 
lnf°  must  be  a  linear  combination  of  the  four  conserved 
properties  'P1,  i  =  1  -*  4,  such  that 


F.  5 


for  m  =  l->3;  i,j,k,l  =  1  -»  5  and  where 


ft '  and  is ‘ 

a  five  vector. 


are  scalars, 
Def ine : 


is  a  three  vector,  and 
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so  that  equation  F.5  becomes 


F 


Equation  F.8  may  be  recast  as: 


4,6  o 

where  ft,  ••  ft*  =  0 .  Thus : 


F 


F. 


Def ine : 


x.»  p.-ip'-wr 


By  Enskog's  method,  the  average  value  of  momentum  is  defined  as 


F. 


Substituting  equations  F.10  and  F.ll  into  equation  F.12  yields 


Therefore : 


Now: 


From  Appendix  C: 
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Finally,  the  temperature  shall  be  defined  as: 

This  is  the  usual  kinetic  theory  definition  of  temperature 
Thus : 


F*  31 

5nWT.  _L  (  R*  Jp  F.  32 
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which  by  Appendix  C  may  be  recast  as: 


Thus,  when  the 
above , 
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then 


Thus : 


and  f°  is: 
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values  of  f°  and  <x/°*  are  substituted  into  the 
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The  desired  equilibrium  properties  are  calculated  from  the 
partition  function  of  the  gas.  Partition  functions  exist, 
however,  only  when  energy  is  quantized.  They  are,  therefore, 
outside  the  slope  of  classical  mechanics.  This  is  a  problem  with 
classical  mechanics  as  a  whole  and  is  not  peculiar  to  the  present 
analysis.  Therefore,  the  works  of  others  shall  be  followed  by 
assuming  a  classical  partition  function  exists  and  is  defined  by: 


Jillf 


For  the  present  case  of  a  diatomic  gas,  g  is  five.  Thus: 


-R ll2tV\ 


where  V  is  the  five-space  volume.  Using  the  results  of  Appendix 
F,  the  classical  partition  function  is: 


The  specific  heat  at  constant  volume  and  the  specific  heat 
at  constant  pressure  are  derived  from: 
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where  U  is  the  internal  energy  of  the  gas.  The  internal  energy 
may  be  expressed  in  terms  of  the  classical  partition  function  as 


U--  NkT V  J2\ 
i  v  «rr  /v 

G.  7 

Thus  : 

a-  mill'  .  X t3/*  =  £  NbT 

-j-'/j  2  2 

G .  8 

which  yields 

--  f  tik 

G .  9 

and 

G.  10 

for  the  equilibrium  specific  heats.  The  ratio  of  specific  heats 
is : 


*  =  SL  ,  I.I00 

Cy 


G.ll 


Expressed  in  terms  of  the  classical  partition  function,  the 
specific  entropy  is: 


s’Ntl'  ‘Mil'll] 


G .  12 


Thus,  substituting  for  Z  yields: 


S  -  Nk 


'•‘•Iftff-'fl-l] 


G.  13 


The  volume  in  five  space  may  be  expressed  as  the  product  of  two 
volumes 
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where  V3  is  the  typical  three-space  volume  and  V2  is  the  two- 

space  'volume'  of  the  angular  coordinates  9  and  0.  Thus, 
equation  G.13  is: 


The  volume,  V2 ,  is  simply  4ir  -  the  solid  angle  of  a  sphere. 


Thus : 


As  known  from  quantum  mechanics,  G.16  is  correct  for  a  gas 
composed  of  heteronuclear ,  rigid-rotor  molecules.  Quantum 
mechanics  predicts  that  in  general 


where  o  is  the  symmetry  number  such  that: 


a  =  1  for  heteronuclear  molecule 


a  =  2  for  horaonuclear  molecules. 


The  symmetry  number  arises  out  of  the  requirement  that  the  wave 
function  of  a  homonuclear  molecule  possesses  a  certain  symmetry 
with  respect  to  the  interchange  of  the  nuclei  of  the  molecule. 
Classically,  the  symmetry  number  ensures  that  one  orientation 
identical  to  another  orientation  does  not  contribute  to  the 
entropy.  A  rigid  rotor,  for  example,  has  identical  orientations 


Since  f  can  be  expressed  as  a  function  of  q^ ,  and  t 
Boltzmann-type  equation  may  be  recast  as: 
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The  Enskog  expansion  of  <JD'  yields 


JD$ 


a  ii*  -  a  ii‘  ^ . 


As  outlined  in  the  text: 
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where  p  is  the  thermodynamic  pressure  defined  as  p  =  nkT. 
Substituting  the  above  into  H.2  gives 
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Expanding  the  first  right-hand  term: 
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Dti  J  Dt  Qfc 


and  substituting  for  from  the  development  of  the  text  yields 


Substituting  for  2«i  from  the  text  gives 
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To  the  first  approximation  then,  the  variation  of  temperature 
during  the  motion  of  the  gas  follows  the  adiabatic  law: 
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Equation  H.7  becomes  then: 
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Substituting  for  yields: 
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Dt 

The  sum  of  the  first  and  last  terms  of  equation  H.3  is: 
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The  bracketted  term  is  a  non-divergent  tensor  which  shall  be 

D*n 

designated  as  'i.u/B,j  .  Thus,  the  above  equation  becomes: 


The  second  and  third  terms  of  equation  H.3  are: 
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which  by  the  definition  of  the  thermodynamic  pressure  is 


The  bracketted  term  is: 
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which  upon  substitution  into  equation  H.20  yields 


'‘'"'I’ 


M  JT.2..M] 


Equation  H.3  is  thus: 


Dl'-V 


BiL Usitf  Jji  J  BlVTL  Bi 


■ffl 


H.  25 


By  Enskog's  solution  method: 
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Thus : 
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Following  Chapman  and  Cowling1,  the  peculiar  momentum  shall  be 
non-dimensionalized  to  the  reduced  momentum.  The  reduced 
momentum  is  defined  as: 


H.  27 


Defined  in  terms  of  the  reduced  momentum,  the  equilibrium 
distribution  function  is: 
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while  equation  H.26  is: 
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Now: 
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for  1  =  1,4,5 

when  1  =  1,4,5  because  the  above  is  evaluated  at  a  point  in 

space.  The  second  integral  is  zero  for  all  'F1  because  the 

integrand  is  either  an  odd  function  of  ?  or  it  is  a  non-divergent 

tensor.  For  f1  and  ¥4  the  first  integral  vanishes  because  the 

integrand  is  an  odd  function  of  ?.  For  !F5 ,  the  first  integral 
may  be  recast  as: 
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where  A  is  a  constant.  The  solution  of  H.35  is: 
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Thus,  the  conditions  of  solvability  for  the  Enskog  method  are 

satisfied.  Therefore,  since  the  general  solution  is  the  sum  of 

the  particular  solution  and  the  homogeneous  solution: 

3 
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2,1 

where  3 

is  a 

three  vector  and  ^ 

is  a  five-vector  where 
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By  identity,  we 

may  define 
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Thus , 
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Since  equations  H.31  and  H.32  are  functions  of  n,  and  T 
only,  then  by  arguments  similar  to  those  of  Chapman  and  Cowling1: 
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where  J((?)  and  K(ff)  are  scalar  functions  of  (? ,  n,  and  T.  To 
satisfy  the  requirement  that  the  average  molecular  quantities  be 
solely  related  to  f  ,  we  must  set: 


jtTT 


jp  =  C> 


H.  43 


for  i  =  1,4,5 


Thus,  neglecting  vanishing  integrals  yields 

i  =  1 
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i  =  4 


JneV-'pi^il  js=c 
28: 


H.  45 


i  =  5 


H.  46 


Equations  H.44  and  H.45  imply  that  ^  and  (3'  are  both  zero. 

1  i  *' 

Equation  H.46  implies  that  (*4  and  ^  are  equal  and  opposite  in 
magnitude  for  j  =  1,2,3.  For  j  =  4,5  ,  and  therefore 


4' 


(3j  p*  w  i*  l-*s ;  i*  i.S 


is  proportional  to  T  and  can  be  incorporated  into  the  first  term 
of  the  integrand  of  H.46.  Thus,  we  may  write: 
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Thermal  Flux  Vector; 


By  definition: 


1.1 


Substituting  for  f*  : 


1.2 


1.3 


because  the  integrals  of  odd  functions  of  p  vanish.  Substituting 
for  yields 


1.4 


Equation  1.4  may  be  simplified  using  the  results  of  Appendix  D: 
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By  definition: 


by  neglecting  vanishing  integrals.  Substituting  for  F 


1. 10 


1. 11 


I. 11  is  equivalent  to: 
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From  the  results  of  Appendix  C,  equation  1.12  is: 


which  by  equation  H.42  is 


*  mf 


In  the  notation  of  Chapman  and  Cowling  : 


where : 


Since  the  stress  tensor  has  been  defined  as 


P.^,  equation 
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we  have  for  the  second  approximation: 
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We  have  for  A: 


j.i 

5* i  J  1 

A  linear  combination  of  equations  J.I  and  H.47  yields: 

which  by  equation  H.31  is  equivalent  to: 

01'  J. 

Summarizing  for  a  moment,  the  coefficients  of  heat  conductivity 
and  viscosity  for  a  gas  composed  of  rigid-rotors  are: 

The  application  of  Sonine  polynomials  to  the  above  integrals  is 
very  similar  to  their  application  in  the  case  of  a  monatomic  gas 
Thus,  only  the  minor  arithmetic  changes  are  noted  below.  The 

derivation  again  closely  follows  that  of  Chapman  and  Cowling1. 


From  equation  J.2,  aQ  =  0  and 


where  the  subscripts  on 
since : 


^  have  been  deleted  for  notational 
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... 


we  have: 


The  results  of  Appendix  C  lead  to: 
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which  leads  to 
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Thus : 


Lastly,  the  Prandtl  number  is  defined  as: 


rf  a 


For  a  diatomic  gas  molecule,  Cp  =  3.5k.  Substituting  this  into 
the  above  along  with  equations  J.ll  and  J.19  yields 

P  -  —  j.2 

'  5U 

where  the  higher  order  terms  of  the  Sonine  polynomials  have  been 
neglected. 


REFERENCES 


1.  Chapman,  S.  C.,  and  Cowling,  T.  G. ,  The  Mathematical  Theory 
of  Non-Uniform  Gases,  Cambridge  University  Press,  Cambridge, 

1952. 

2.  Hirschf elder ,  J.  0.,  Curtiss,  F.  C. ,  and  Bird,  R.  B. , 

Molecular  Theory  of  Gases  and  Liquids.  John  Wiley  and  Sons.,  New 
York,  1954. 

3.  McQuarrie,  D.  A.,  Statistical  They modvnamics .  Harper  &  Row, 
New  York,  1973. 

4 .  Bo 1 2 ,  R .  G . ,  and  Tuve ,  G .  L . ,  Handbook  of  Tables  for  Applied 
Engineering  Science,  CRC  Press  Inc.,  Boca  Raton,  Florida,  1981. 

5.  Ishida,  Y. ,  'Kinetic  Theory  of  Rigid  Molecules',  Phys  Rev,  Vol 
10  No  4,  pp  305-321,  1917. 

6.  Pidduck,  F.  B. ,  'The  Kinetic  Theory  of  a  Special  Type  of 
Rigid  Molecule',  Proc  Royal  Society  of  London,  Vol  A101,  p.  101, 
1922. 

7.  Curtiss,  C.  F. ,  'Kinetic  Theory  of  Nonspherical  Molecules', 
Journal  of  Chemical  Physics,  Vol  24  No  2,  pp  225-241,  1956. 

8.  Kirkwood,  J.  G. ,  'The  Statistical  Mechanical  Theory  of 
Transport  Processes:  I.  General  Theory',  Journal  of  Chemical 
Physics,  Vol  14  No  3,  pp  180-201,  1946. 

9.  Kirkwood,  J.  G. ,  'The  Statistical  Mechanical  Theory  of 
Transport  Processes:  II  Transport  in  Gases',  Journal  of  Chemical 
Physics,  Vol  15  No  1,  pp  72-76,  1947. 


103 


U. 5. Government  Printing  Office:  2  986  18  06] '60 


